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Abstract 

We propose a theory of semiclassical mechanics in phase space based 
on the notion of quantized symplectic area. The definition of symplectic 
area makes use of a deep topological property of symplectic mappings, 
known as the "principle of the symplectic camel" which places stringent 
conditions on the global geometry of Hamiltonian mechanics. Following 
this principle, symplectic mappings -and hence Hamiltonian flows- are 
much more rigid than Liouville's theorem suggests. The dynamical objects 
of our semiclassical theory are "waveforms" , whose definition requires the 
notion of square root of de Rham forms. The arguments of these square 
roots are calculated by using the properties of a generalized Maslov index. 
The motion of waveforms is determined by Hamiltonian mechanics, and 
the local expressions of these moving waveforms on configuration space 
are the usual approximate solutions of WKB-Maslov theory. 
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1 Introduction 

Non-relativistic physics is governed by two Sciences with distinct domains of 
apphcabihty: classical mechanics (CM), and quantum mechanics {QM). The 
paradigm of CM is Newton's second law 

while that of QM is Schrodinger's equation 

dm 

ih— = H^. 

dt 

Both equations describe motions (of particles in the first case, and of waves 
in the second) in configuration space R". However, CM and QM differ pro- 
foundly, both physically and mathematically. They differ physically, because 
QM renounces to the idea of material systems with sharply defined positions 
and momenta, and incorporates instead complex probability amplitudes in its 
dynamics. They differ mathematically, because while Newton's second law can 
immediately be interpreted in terms of phase space variables using the Hamil- 
tonian formalism, there is no simple and obvious way to define "phase space 
wavefunctions" . On the other hand, one of the most useful manifestations of 
QM, both in physics and chemistry, is semiclassical mechanics (SM), which 
applies when the scale relative to h of certain parameters (e.g. position, time, 
or mass) in a system is large. Systems to which SM applies exhibit a behavior 
which is both classical and quantal: while certain quantities (for instance energy 
or angular momentum) remain quantized, the motion of the system is governed 
by CM. (SM is sometimes described as a way of doing a simplified path-integral 
formalism with a focus on "classical paths".) 
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The aim of this article is to present a unifying and mathematicahy rigorous 
theory of semiclassical mechanics in phase space based on a deep and striking 
topological property of Hamiltonian flows, the non-squeezing theorem. This 
theorem -also known as the "principle of the symplectic camel"- says that no 
Hamiltonian flow will ever be able to squeeze a phase space ball into a phase 
space cylinder of smaller radius based on a plane of conjugate variables Xj,pj. 
We will use the principle of the symplectic camel, together with a simple physical 
postulate related to the quantization of action, to quantize phase space in such a 
way that we recover the usual semiclassical energy levels for integrable systems 
by a purely topological argument, without any reference, whatsoever, to the 
WKB method or to approximate wavefunctions constructed by other methods. 

Our paper consists of two parts, which can be read independently: 

• In the first part (Sections ^ and|3|) we begin by reviewing the "principle of 
the symplectic camel" (which seems to be little known by physicists) . We 
define the related notion of symplectic area, which we then use to quantize 
energy shells by an appropriate physical postulate on the periodic orbits 
they carry. This postulate remarkably leads to the correct ground energy 
levels for the anisotropic harmonic oscillator in arbitrary dimensions (it 
can also be used to derive a classical form of Heisenberg's inequalities as 
we have shown in We show that our postulate leads, by a topological 

argument, to the usual Keller-Maslov quantization condition 



in the integrable case (7 a loop on the "invariant torus" ) . Our quantization 
procedure is actually much more general than those find in the literature 
(it quantizes periodic orbits and energy shells) , and could thus be applied 
with profit to systems exhibiting chaotic behavior. 

• The second part (Sections IV to VII) begins by the study of a simple 
example, the one-dimensional harmonic oscillator. It contains in embry- 
onic form the whole theory which is being further developed in Sections 
V to VII. We then proceed to survey the notion of phase of a Lagrangian 
submanifold, as defined by Leray and the cohomological theory of the 
Maslov index which we have developed in |l|, . We are thereafter 
able to define our "waveforms" : they are phase objects whose phase is ex- 
pressed in terms of the universal covering of the Lagrangian submanifold, 
and whose amplitude is the square root of an arbitrary "twisted" (or de 
Rham) form. Our study of the Maslov index will allow us to assign the 
proper argument to these square roots. Our constructions apply whether 
the underlying manifold is oriented or not (in contrast with other quan- 
tization theories where orientability is a sine qua non requirement, as for 
instance in [^). Finally, we show that the local expressions on configu- 
ration space of our waveforms, whose motion is Hamiltonian, are just the 
usual WKB wavefunctions. 




pdx 



—771(7) is an integer 
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Notations and terminology. 



The letter z denotes the generic point {x,p) of the phase space M^" — M" x 
Kp. We equip K^" with the standard syniplectic form = d{pdx): 

n 

Q — dp A dx = dpj A dxj . 
i=i 

We will denote by A(n) (resp. Sp{n)) the Lagrangian Grassmannian (resp. the 
symplectic group) of the symplectic space (K.^", O): ^ G A(n) if and only if £ is 
a n-dimensional linear subspace of R^" having the property that il{z, z') = 
for all z, z' . The symplectic group Sp(n) consists of all automorphisms s of 
R^" such that s*fi = ^2, that is ri(sz, sz') = ri(z, z') for all z, z' . The universal 
coverings of A(n) and Spin) will be denoted by Aoo('t-) and Spoa[n\ 

By definition, a symplectic transformation (or: canonical transformation) is 
a diffeomorphism of phase space whose Jacobian matrix belongs to Sp{n) at 
every point at which it is defined. Also recall that a Lagrangian submanifold 
of phase space R" x R^ is a n-dimensional submanifold V C R" x R^ whose 
tangent spaces are all Lagrangian planes. Equivalently, iyfl = 0, Ly being the 
inclusion operator F C R" x R". 

We will also use in Section O some elementary notations from singular (co-) 
chain theory. Let X be a non-empty set, (G, +) an Abelian group and p an 
integer > 0. A G-valued p-cochain on X is a mapping c : A'^'+^ G. The 
coboundary of c is the {p + l)-cochain dc defined by 

P+2 

dc{xo,...,Xp+i) = ^{-ly c{xq, ...,Xj, ...,Xp+2) 

where the cap " deletes the term it covers. If dc — 0^ c is called a p-cocycle; 
if c = 96 for some {p — l)-cochain, it is called a coboundary. We have = 0, 
hence a coboundary is a cocycle. 

2 Symplectic Camel Quantization 

Let B{R) be a closed ball in phase space with radius R: 

B{R) = {{x,p) : |x- + b-poP < R^] 

and Zj{r) a cylinder with radius r: 

Zj(r) = {ix,p) : (xj - XQ^j f + {pj - Pa,jf < r^} 

(1 < J < n) based on the Xj,pj plane (we will call hereafter the Zj{r) symplectic 
cylinders) . Gromov pl| proved in the mid 1980's that there cannot exist a sym- 
plectic transformation sending B{R) inside Zj{r) unless R < r. In particular. 
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a Hamiltonian flow can never squeeze a phase space ball inside a symplectic 
cylinder with smaller radius. Gromov's theorem is equivalent statement of the 
principle of the symplectic camel : 

Proposition 1 Let Pij : x Mp — > R^.^ x Mp^. be the projection operator. 
For every symplectic transformation f , we have: 

AreaPr J f{B{R)) > irR^. (1) 

(See the proof in where we used (|l|) to derive a classical form of Heisen- 
berg's uncertainty relations.) 

This result is of course striking, because it seems to contradict the common 
conception of Liouville's theorem, which is that under a Hamiltonian flow a 
volume in phase space can be made as thin as one likes (cf. Gibbs who calls 
this the "principle of extension in phase" ; also see the discussion of Liouville's 
theorem in Penrose [0). However, what is overseen is that the proof of Liou- 
ville's theorem only uses the fact that Hamiltonian flows are divergence free. 
In fact, Hamiltonian flows consist of symplectic transformations, and this is a 
much stronger property than being just volume preserving as soon as n > 1. 
For instance, the statement 

" / is a volume-preserving transformation of phase space^^ 

is equivalent to saying that if {x,p) = f{x',p') then the Jacobian matrix 

( dx dx \ 
dx' dp' / 

has determinant equal to one, while the statement 

f is a symplectic transformation of phase space^^ 
means that the entries of f'{z) satisfy the much more stringent conditions 

iMf§§^^ (^j'^^ are symmetric, 

f d^\'^ dp_ _ f dp_\ dx _ J 

\ dx') dp' \dx' ) dp' — ^n-Kn- 

No "easy" proofs of Gromov's theorem are known. In Gromov's original 
paper and in Hofer-Zehnder the reader will find proofs making use of the 
theory of pseudo-holomorphic curves. Viterbo gives in | |3^ a very interesting 
alternative proof using the notion of generating function. 

2.1 Symplectic Area and Periodic Orbits 

Let I? be a subset of R" x R^ . We will call symplectic radius of T) the supremum 
^max of all i? > such that we can send the phase space ball B{R) inside V 
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using a symplectic transformation. We will call symplectic area of 2?, and denote 
by A{'D) the number ttR^^^: 



A{V) 



sup {7ri?2 : f{B{R)) C V} . 



f symplectic 



{A{'D) is also sometimes called the symplectic capacity of D; see e.g. p3[). It 
immediately follows from the definition of .4(1?) that — AiT)) for every 

symplectic transformation / of phase space: symplectic area is thus a symplectic 
invariant. 

Remark. The notion of symplectic area was first introduced by Ekeland 
and Hofer in Q; there are other non-equivalent definitions of symplectic ar- 
eas/capacities (see e.g. [p3|). 

The principle of the symplectic camel can obviously be restated as 



showing that subsets of phase space with very different shapes and volumes can 
have the same symplectic area. Let for instance 



be a phase space ellipsoid; we assume that Ri < ■■■ < Rn- (The equation of every 
ellipsoid in phase space can be put in the form above by a suitable symplectic 
change of coordinates.) It follows from property (^) that the symplectic area of 
this ellipsoid is 



The symplectic area of a set has -as the terminology suggests- the dimension 
of an area. In the case n — 1 the symplectic area is in fact just the usual area: 



Notice that the symplectic area of a ball B{R), or of a symplectic cylinder, is 
independent of the dimension of the ambient phase space, as it always is ttR'^ . 
The symplectic area and the volume of a ball B{R) in M" x are related by 
the formula 




(2) 




A{E) ^ttRI. 



(3) 




(4) 



Yo\B{R) = [A[B{R))] 



n 



(5) 



since B{R) has volume 7r"i?^"/n!. 



6 



2.2 Symplectic Area and Periodic Orbits 

Consider a bounded domain T) with boundary 7 in the phase plane 
Obviously formula (0) can be written, using Stoke's theorem as 



pdx 



1 



A{V)- 
showing that 

symplectic area — action" 

in the case n = 1. It turns out -and this is another striking feature of 
the principle of the symplectic camel- that this relation holds in any dimen- 
sion. In fact, symplectic area is related to the action of periodic orbits of 
Hamiltonian systems. Let us begin by some general considerations. Consider 
an infinitely differentiable function ("Hamiltonian") H : R" x — > R and 
Xh = {VpH, —VxH) the associated Hamilton vector field. The periodic orbits 
of Xh are defined as follows: let {ft) be the flow of Xh and assume that there 
exists z and T > such that /t(z) = ft+riz)- Then 7(i) = /t(z), < i < T is a 
periodic orbit through z. The action of the periodic orbit 7 is then the integral 

pdx = / p{t)dx{t) 
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where {x{t),p{t)) = ft{z). 

By definition an "energy shell" E of is. a non-empty regular level set of 
the H: 

E = {{x,p) e R^ X R;' : H{x.,p) = E] . 

Any smooth hypcrsurface of phase space can of course be viewed as the energy 
shell of some Hamiltonian function H: it suffices to choose for H any 
function on R" x R^ keeping the constant value E va a, tubular neighborhood 
of S. By definition, a periodic orbit of E is then a periodic orbit of the flow 
determined by H , and lying on the energy shell E. Of course, for this deflnition 
to make sense, we have to show that these periodic orbits are independent of 
the choice of the Hamiltonian having E for energy shell. This follows from the 
following well-known result: 

Lemma 2 let H and K be two Hamiltonians, and suppose that there exist two 
constants h and k such that 

E = {z : H{z) = h} = {z: K{z) = k} (6) 

with ^ and 'VzK 7^ on E. Then the Hamiltonian vector fields Xh and 
Xk have the same periodic orbits on E. 
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Proof. It suffices to show that Xh and Xk have the same integral curves, up 
to a reparametrization. Since VzH(z) ^ and Vz-ftr(z) ^ are both normal 
to E at 2, there exists a function a 7^ on E such that Xk = aXn on E. Let 
now (/f) and [jjt) be the flows of H and respectively, and define a function 
t = ^(z, s), s € K as being the solution of the ordinary differential problem 

^=a(/,(z)) , i(z,0) = 
as 

(where z is being viewed as a parameter). We claim that 

as{z)^h{z) for zeE. (7) 

In fact, by the chain rule 

^h[z) = tM^)t = XH{ft{z))a{h{z)) 
as at as 

that is, since Xk = aXn- 

±f,[z)^XK{ft{z)) 

which shows that the mapping s 1 — > ft(z.s){z) is a solution of the differential 
equation z = X}j{z) passing through z at time s = t{z, 0) = 0. By the unique- 
ness theorem on solutions of systems of differential equations, this mapping 
must be identical to the mapping s 1 — > gs{z)] hence (^. Both Hamiltonians H 
and K thus have the same periodic orbits. ■ 

The general problem of the existence of periodic orbits on a given energy 
shell E is a very difficult one, which has not yet been completely solved. We 
have however the following partial result (see |Q and the references therein): 

Proposition 3 // the hypersurface E is the boundary of a compact star-shaped 
submanifold of phase space, then it carries at least one periodic orbit. 

(Recall that a submanifold M of an Euclidean space is called star-shaped if 
there exists a point z G M such that the line segment joining z to any other 
point z' G M lies inside M.) In particular, the boundary of every closed convex 
submanifold thus carries a periodic orbit. 

The essential relation between the action of periodic orbits and symplectic 
area is given by the following theorem: 

Theorem 4 Let M be a compact star-shaped submanifold in phase space. Then: 
(1) Every periodic orbit ^ onY, = DM is such that 



pdx 



> AiM) (8) 



and: (2) There exists at least one periodic orbit 7jjji,j whose action is the sym- 
plectic area of M : 



pdx 

Tmin 



= A{M). 
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(See again Hofcr-Zehnder's treatise |Q for a proof.) 

Remark. We conjecture that the property of the symplectic camel is the 
key to a better understanding of not only quantum mechanics, but also of clas- 
sical phenomena. Consider, for example, adiabaticity. While it is rather well 
understood in one dimension (cf. "Einstein's pendulum"), one must take the 
usual physical statements and "proofs" of adiabatic invariance in higher dimen- 
sions with more than a critical eye. The existence of the symplectic invariant 
A{M) shows that symplectic geometry is, in a sense, a two-dimensional geom- 
etry "projected" in higher dimensions. Perhaps a general adiabatic principle 
could be derived from Proposition |l| by showing that adiabatic invariance in 
the phase plane is sufficient for deducing more general results. We also conjec- 
ture that the principle of the symplectic camel might play a fundamental role 
in thermodynamics and statistical physics (e.g., Bose-Einstein and Fermi-Dirac 
statistics). Viterbo |Q has given other interesting physical interpretations of 
the principle of the symplectic camel. 



3 Phase Space Quantization 

The property of the symplectic camel discussed above can be used to quantize 
phase space in a very simple way. We will, in particular, recover the correct 
ground energy level for the n-dimensional anisotropic oscillator. 

3.1 A Physical Postulate 

We now make the following postulate of physical nature: 

Minimum symplectic area postulate: The only physically admissible 
periodic orbits are those which lie on hypersurfaces E which are boundaries of 
convex subsets M of phase space with symplectic area at least ^h. Moreover, if 
A{A-I) — \h then E effectively carries a physically admissible minimal periodic 
orbit 7„,i„. 

Notice that in view of Theorem ^ such a minimal periodic orbit satisfies 

pdx^\h (9) 

hence our postulate implies quantization of action. It actually implies much 
more, as we are going to see: because of the principle of the symplectic camel, 
it is not a mere restatement of (^). As we remarked in the Introduction to 
this paper, our Postulate quantizes energy shells and periodic orbits, and thus 
applies beyond integrable systems. 

We begin by giving an immediate striking application, by showing that the 
minimum symplectic area postulate leads to the correct energy levels of the 
anisotropic multi-dimensional harmonic oscillator. 
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Proposition 5 Consider the n- dimensional harmonic oscillator with Hamilto- 
nian: 

n ^ 

^ = E2^(p|+-.-H) (10) 

The minimum symplectic area postulate implies that the ground energy level of 
that Hamiltonian is 

n 

Eo = Y.\t^j. (11) 

Proof. Let L be the diagonal matrix with diagonal entries {mjLUj)~^^^ . The 
symplectic change of variables {x,p) i — > {Lx,L~^p) changes H into 

n 

The change of variables preserving both action integrals and symplectic areas, 
it is sufficient to prove the theorem for H' . Each orbit 




xi = x[ cosuit + p[ sinwif , pi = x[ sinivit — p\ coswif 



Xn = x'^ coswni sinw^i , p„ = a;^ sina;„i - p'^ cosw„t 



lies, not only on the ellipsoid which is the energy shell of the Hamiltonian W , 
but also on each of the symplectic cylinders 

Z^(Rj) = {{x,p):x']-rp] = Rfi 

with i?? = x'J + p'J and 1 < j < n. These cylinders carry periodic orbits, and 
their symplectic areas must thus satisfy the conditions 

in view of our postulate. If 7mi„ is a minimal periodic orbit, it will thus satisfy 

n n 

which is the result predicted by standard quantum mechanics. 
3.2 Quantization of Integrable Systems 

Let us next consider a completely integrable system with Hamiltonian H. There 
are thus n independent constants of the motion Fi = iJ, F2, F„ in involution: 
{Fj,Fk} = for 1 < j,k < n. Given an energy shell S of i?, through every 
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point zq — {xo,Pa) of S passes a Lagrangian submanifold V carrying the orbits 
passing through zo- Moreover, when V is connected (which we assume) there 
exists a symplectic transformation 

f -.V — > {S^)'' X K"-'^ (12) 

where {Sj}'^ is the product of k unit circles Sj, each lying in some coordinate 
plane Xj,pj (/ can be constructed using "action-angle variables", see e.g. [Q, 
p^). The minimum symplectic area postulate imposes a condition on the energy 
shells of any Hamiltonian. That condition is that there should be no periodic 
orbits with action less than i/i, and that there should exist "minimal periodic 
orbits" having precisely as action. In fact, we have the following result which 
ties the minimum symplectic area/action principle to the Maslov index of loops: 

Theorem 6 Let V be a Lagrangian submanifold associated to a Liouville in- 
tegrable Hamiltonian H and carrying minimal action periodic orbits. Then we 
have 

^J^vdx-\mi,)^, (13) 

for every loop on V . 

Proof. Since the actions of loops are symplectic invariants, we can use the 
symplectomorphism (^2|) to reduce the proof to the case V = (S'^)'^ x M"^'''. 
Since the first homotopy group of V is 

TTliiS^)'' X W'-^) EE TTiiS^)'' = (Z^+) 

it follows that every loop in V is homotopic to a loop of the type: 

7(0 = (7i(i), -,7^(0, 0,...,0) , 0<t<T 

where 7^ are loops on S^: 7^(0) — jj{T). On the other hand, every loop on 
is homotopic to a loop ej{t) ~ {cos ojjt, sin ujjt), < t < Tj so that there must 
exist positive integers Hj (1 < j < n) such that fJ-iTi — ■ ■ ■ = fJ-i^T^ = T. We 
can thus identify 7^ with iJ.jSj, the loop ej described "/i^ times": 

IJ.jej{t) — {cos Ldjt^siwLOjt) <t <T 

and it follows that any loop in V — {S^) ^' X R"-'= is homotopic to a loop 
7 — ^i^ei + • ■ • + /ife£fe- We thus have 

k 



(j) pdx — /ij (j) pj dxj 

J-l j^i Jej 



and using the same argument as that leading to the proof of formula ( |ll| ) in 
Proposition we must have 



<i) Pjdxj — \h (1 < J < k) 
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and hence 




Now, the Maslov index of such a loop 7 in by definition 

fc 

(see formula (|5|) in Example |l^, Section ^) hence the Keller-Maslov condition 
(©■ ■ 

This result motivates the following definition: 
Definition 7 A Lagrangian submanifold V is said to be quantized if 

pdx — —m{'-f) is an integer (14) 

for every loop ^inV. 

This definition is of course nothing else than the usual Maslov-Keller quan- 
tization condition ||2^, originating historically from WKB theory. We 
arrived to it by purely topological considerations. 




4 Waveforms on the Circle 

We consider in this Section the one dimensional oscillator with Hamiltonian 
function 

The flow determined by Hamilton's equations for H consisting of the rotations 

_ / cost sint \ 
\^ — sint cost j 

the phase-space trajectories are thus the circles 5*^ = {|z| = r}. These circles 
carry a natural length element denoted, with the usual abuse of notation, ds = 
rdO, where 9 is the polar angle. 



4.1 Position of the Problem 

One wants to define on S}. objects whose vocation is to play the role of waveforms 
in phase space, in the sense that their local expressions are, at best the "true" 
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wavefunction, at worst their semiclassical approximation (i.e. approximations 
for "small fi," ) . One looks for an expression of the type 

^{z) = eTi'^^''>a{z)Vds (15) 

where the phase ip and the amplitude a are real functions, and ^/ds is supposed 
to have some well defined meaning . Unfortunately, one immediately encounters 
two difficulties when one tries to define (p and \/ds. First of all, if one wants 
the theory to be consistent with semiclassical mechanics, one must require that 
the differential of the phase be the action form: 

dip = pdx = -r^ain^ ede. (16) 

Unfortunately there exists no such function ip because the 1-form pdx is not 
exact on . We can however define a function (p satisfying (^6|) on the universal 
covering tt : M. ^ oi S^. That covering is defined by 7r(6') = r{cos9,sm9) 
and one immediately checks that 

Vi&) = Y (^^^ 6' cos 61 - 0) (17) 

satisfies (|l6|). We are thus led to consider ^'(z) as being an expression of the 
type 

^{9) = eTi'f'^'^'^ a{9)V7d0 

where one allows 9 to take any real value, which amounts to define the candidate 
for being a phase space wavefunction on the universal covering of the circle. 
However, there is a second, more serious obstruction because one does not see 
how to define unambiguously the square root \/ds — \frdB. The simplest way 
out of this difficulty is to decide that one should only consider the (for instance, 
positive) square root of the density \ds\, that is that we take 

^(0) = et'^(^)a(0)^7d^ (18) 

which indeed has a well defined meaning. However, there is a serious rub with 
that choice because it leads to the wrong energy levels: since we are actually 
interested in a single- valued function on S^, we have to impose the condition 

^'(61 + 27r) = *(6') (19) 

to the expression (|l8|), which is equivalent to the condition 

ip{e + 2n) ^ ip{e) ~ 2NTTh 

for some integer N. By definition of (|l^ this is in turn equivalent to — 2Nh, 
which leads to the energy levels En = Nh, instead of the physically correct 
En = {N+^)h. 
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4.2 The Need for De Rham Forms 



The way out of these difficulties, and which leads to the correct quantization 
conditions, is the inclusion in the theory of de Rham's Q "forms of odd kind" 
(also called "twisted" or "pseudo" -forms in the literature) are to ordinary dif- 
ferential forms what "pseudo- vectors" are to ordinary vectors. By this we mean 
that the local expressions of the de Rham forms depend on the orientation of the 
local charts used to define them (rigorously speaking, the de Rham forms are 
just ordinary forms, but defined on the oriented double cover of the manifold). 
In the case of the harmonic oscillator, this leads to the following constructions. 
Consider the atlas of consisting of the four half-circles 

S'r,up = : = r, Imz > 0} , 5i,_„ = {z : \z\ = r, Imz < 0} 
Sr,ieft = U ■ \z\ = r,Re2; < 0} , Sl„g^t = {z : \z\ = r,Rez > 0} 

together with the projections /„p,/dou.« : {x,p) — > x and fieft, fright ■ {x,p) — >■ 
p. The atlas thus defined is not oriented; for example, the transition function 
on Si jgjj n 5*^ has negative sign. The local expressions of dO in the charts 
defined above are, respectively 



(d0)„p=£(r2-a;2)-l/2rf2. ^ (^0)^^^^^^ ^ ^(^2 _ ^2)-l/2 

{de)i,ft = e(r2 - p2)-i/2 ^ [deUght = e{r^ - p")-^'^ dp 



(20) 



where e = ±1 is the orientation induced from the x, p axes on 5*^^, etc. by 
the diffeomorphisms /„p, etc. Thus, if the axes come equipped with their usual 
orientations, then e = —1 for {d9)up and {dO)ieft and -1-1 for {dO)down and 
{d6)right, and a change of orientation has the effect of reversing the sign of e. 
The formulas (pO|), which are characteristic for de Rham forms, suggest that we 
define the "argument" of d9 by 

r m{e)Tr in 5^ \ {±r} 
argde^l (21) 
[ {m{e) + l)7r in S^ \ {±ir} 

where the integer m{9) is defined by 

m{9) = [e/n] + 1 (22) 

the square brackets meaning "integer part of" . Notice that a change of orienta- 
tion of the frame x,p amounts replacing m(9) by m(9 + 7r) — m{9) + l. Formulas 
(|l]) and (|2^) allow us to define the square root of ds ~ rd9 in each of the sets 
sf\ {±r} and Si \ {±ir}. In fact. 



ds = i'^^^'^ ./\FdB\ in S';t\{±r} 



jm(e)+iy|7^ in Sl\{±ir} 

(Notice that both expressions do not coincide on the overlaps). We are thus 
lead to give the following definition of '^{9): it is the phase space object whose 
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expression on \ {±r} is given by 

*o(^) = e-^'^(^)a(0)i"W^7d^ (23) 
and on \ {±ir} by ^'1(6') = i*(6i): 

^i{9) = ei'^('')a(0)i"(«)+i^7d0j (24) 

With that definition the single-valuedness condition ( p^ ) becomes 

^',(^^ + 2^) = , J = 1,2 (25) 

and is equivalent to = {2N + l)h, which yields the true energy levels En = 
{N + ^)fi predicted by quantum mechanics. 

We are going to show that this construction of phase space waveforms can 
be extended to any physical system to which a Lagrangian submanifold can be 
associated. We begin by defining a notion of phase on arbitrary Lagrangian 
manifolds generalizing (p7|). 

5 The Lagrangian Phase 

In the rest of this article the letter V will denote a connected Lagrangian sub- 
manifold. Lagrangian manifolds are associated in a natural way to integrable 
classical physical systems, and to every quantum system: 

Example 8 The integrable systems of classical mechanics: V is then topologi- 
cally an "invariant torus or, more generally a product of k circles and n — k 
lines. 



Example 9 Let ^p{x) = a(a;)e"s^*^^^ where a and $ are defined on some con- 
nected open subset of configuration space. The graph V = {{x,\7 x^{x)} is a 
Lagrangian submanifold. 

5.1 Definition of the phase 

Consider the universal covering tt : V — > V of the Lagrangian submanifold V. 
Since V is simply connected there exists a differentiable mapping ip : V — > M 
such that 

d(f{z) — pdx if 'n:{z) — {x,p). (26) 

We will call, following Leray |^, such a function Lp a phase of V . The phase 
can be explicitly constructed in the following way: choose an "origin" zq G V , 
and identify z ^ V with the homotopy classes (with fixed endpoints) of paths 
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in V originating at zq] the projection Tr{z) is then the endpoint z of an element 
^zoz '^f the homotopy class z. A phase function is then given by the formula 

ip{z) = pdx. (27) 

Clearly the integral only depends on the homotopy class z of J^oz view of 
Stoke's theorem, because il — d{pdx) is zero on V. Also, 

dip{z) = pdx if 7r(z) = {x,p) . (28) 

We observe that the action of the first homotopy group tti{V) = tti{V, zq) on V 
is reflected by the formula 

(p{jz) = (p{z) + pdx (29) 

for all 7 e TTi{V). Thus (p is defined on V if and only if all the periods J^pdx 
of pdx vanish, i.e. if V is contractible. We leave it to the Reader to check that 



formula ( |27| ) leads to the function (17) if we require that ip{0) = 



5.2 The action of Hamiltonian flows on the phase 

Consider a function H = H{x,p,t) defined on some open subset _D x Mt of 
the extended phase space M" x x . We do not assume here that H has 
any particular form (for instance "kinetic energy + potential" ) , but only that 
it is a continuously differentiable function; we also make the simplifying, but 
not essential, assumption that the solutions of the corresponding Hamilton's 
equations 

X^VpH ,p=-W:,H. 

exist for all times, and are uniquely determined by their values at a time t'. 
We denote by {ft,t') the associated time-dependent flow: ft^t' is the symplectic 
transformation that takes a point {x',p') — {x(t'),p{t')) to the point {x,p) = 
{x{t),p{t)). When t' = 0, we write simply ftfi — ft- The time-dependent flow 
satisfies the Chapman-Kolmogorov relation 

ftj'ft',t" = ft,t" (30) 

for all times t, t', t" . 

Suppose that we are given, at some time t', a Lagrangian submanifold Vti, 
and select a base point Zj/on Vt' ■ This allows us to define the phase (p{z,t') of 



Vf by formula (27), with zo replaced by Zf, z being an element of the universal 
covering Vt' ■ The manifold Vt — ft,t'{Vt') is also Lagrangian; defining a base 
point zt G Vt by zt — ft.t'{zt'), we identify the universal covering Vt with Vt', 
defining the projection nt : V ^ Vt' by TTt{z) — z{t) if TTt'{z) = z{t'). Denoting 
by (/5(z, t) the phase of Vf , we have: 
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Proposition 10 The phases ip{z,t) and (p(z,t') are related by the formula: 

l-z(t) 

ip{z,t) ^ ip{z,t') + pdx-Hdt (31) 

Jz{t') 

where the integral is calculated along the trajectory s — > fs,t'{z(t')) (f < s < t) 
leading from z{t') G Vt' to z{t) G Vt- 

Proof. We first note that in view of the relative invariance of the Poincare- 
Cartan form (see pq]), we have 



/•2(t) e rzt 

ip{z,t')+ / pdx-Hdt^ / pdx+ / pdx-Hdt (32) 



where the integral in the left-hand side is calculated along the trajectory s 
fs,t'{zt') {f < s <t), and ft,r{lzoz) is the image in Vt by ft^t' of a path in Vr 
whose homotopy class is z. Denoting by x(z, i) the left hand-side of (32), we 
thus have, for fixed t: 

dx(z,t) = p{t) dx{t) 

so that xiZft) and (p{z,t) differ by a function K{t) only depending on t. Since 
x{z,t') = (p{z,t'), we must have K ~Q. ■ 

5.3 Phase and generating functions 

The notion of phase of a Lagrangian submanifold is related (as is the action 
integral, see [Q) to the notion of generating function. 

Recall (see for instance Arnol'dQ or Goldstein jlo|) that a symplectic trans- 
formation / is free if there exists a function W defined on twice the configuration 
space and such that if {x,p) = f{x',p') then 

pdx — p dx -f dW{x^ x'). (33) 

The function W is then called a, free generating function (or: generating function 
of the second kind) for /. When / is free, the relation {x^p) — f{x',p') uniquely 
determines x' in terms of x. In fact, (p3) being equivalent to 



p^\/^W{x,x') and p' ^ ~\/^>W{x,x') (34) 
we have by the implicit function theorem: 

-|^--^^»- <-) 

Suppose now that iJ is a Hamiltonian function of the type 

n 

H = Y.—{p,-A,{x.t)f + U{x,t). (36) 
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It is then easy to prove (see [|l^) that there exists e > such that for 

0<\t-t'\<e (37) 

the mappings ft^t' are free {{ft,t') is the time-dependent flow determined by H). 
Let now Vf , Vt be as in Proposition Keeping initial position and time x' and 
t' fixed, every point x is thus reached, after time t — t' , by a unique trajectory 
r emanating from x. Suppose now z = {x,p) G Vf. That point is the image by 
ft.t' of a unique point z' = [x' ,p') G Vt> ■ The mapping x i — > x' thus defined is a 
local diffeomorphism, whose inverse we will denote by //^/. Thus, by definition, 
//^/(x') is the unique element of M" such that 

ft,t'{x',p')^{!^Ax').p) ■ (38) 

The action integral is then, by definition, the integral of the Poincare-Cartan 
form along that trajectory; we note that the function 

S{x,x';t,t') ^ / pdx-Hdt (39) 

J x' ,t' 

satisfies Hamilton- Jacobi's equation with initial condition t = t' : 
dS 

— + H{x,V^S,t)^0 , S.,,t'{x,x';t,t) = 0. (40) 

From these considerations we easily get the following consequence of Propo- 
sition |l^: 

Corollary 11 Under the assumptions above on the Jt,t' we have 

ip{z,t) = ip{z,t') + S{x,x';t,t') (41) 

where z has projection Ht'iz) = {x',p') on Vf , and {x,p) — Jt,t'{x' ,p'). The 
local expression 

$(x, t) = $(a;', t') + S{x, x'; t, t') (42) 

ofip{z,t) satisfies Hamilton- Jacobi's equation. 

Formula (fiW is an immediate consequence of (|l|) and (H); formula (H) 
follows from (^0[). (See for a detailed study of the relationship between the 
action integral and free generating functions). 

Remark. When TJ is a quadratic homogeneous polynomial in the Xi^pj, 
Euler's identity for homogeneous functions yields 

H ^]^[x-\/^H +p-\I^p) 

hence, using Hamilton's equations: 

if{z,t) = (p{z,t') -\- ^{p ■ X - p' ■ x') - H{z,t){t - t'). (43) 
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6 The Argument Index 



The construction of an index generalizing the function m(6) ~ [S/tt] + 1 to 
arbitrary Lagrangian manifolds is rather technical and will be done in several 
steps. We have exposed elsewhere (see jl^, |l5|, 16 ) a direct cohomological 



construction of the argument index based on previous work by Leray 
and the author |jl^, ll^. We adopt here a more concrete point of view by 
making use of Souriau's identification of the Lagrangian Grassmannian with the 
manifold of all symmetric unitary matrices (see Souriau's original paper and 
also ||l^, |l^). This approach has the advantage that it allows straightforward 
numerical computations and that it does not require any prior knowledge of 
chain intersection theory. 



6.1 Maslov and Argument Indices for Paths 

We begin by recalling some results from Lagrangian analysis ||l^, |l^, . 
The "Souriau mapping" is the mapping w : A(n) — > U (n) defined by 

w{t) = uu* = u{u^) if ^ = m(M;') (44) 

where u g U(n). This mapping is indeed well defined, because if M(Rp) = u'{W^) 
then u' = uh for some h G 0{n) and hence u'u' = uu* . The mapping w is in 
fact a diffeomorphism, and hence identifies the Lagrangian Grassmannian A(n) 
with the manifold 

W{n) ={w eU{n),w ^ w'^} . (45) 

of all symmetric unitary matrices. The universal covering Aoo{n) of A(n) can 
then identified with the subset 

Woo{n,C) ^ {{w,a):we VF(n, C), det(w) = e'"} (46) 

of U{n,C) X C, the covering mapping being the projection {w,9) i — > w. It 
follows that A(n) = Woo{n)/Z and hence 7ri(A(n)) = (Z, +). The action of 
7ri(A(n)) on Aoc{n) = Woo{n) is given by 

X'' ■ z = {w,a + 2kT:). (47) 

where A is the generator of 7ri(A(n)) whose image in Z is +1. 

Let us write explicitly these identifications in the case n = I. The manifold 
A(l) consists of all straight lines £ through the origin in the phase plane = 
Rj^xRp. We will denote by £{9) the line through the origin whose angle with the 
X-axis is 0+^ (mod tt): £{6) is thus the direction of the tangent to the unit circle 
at the point e^^. Since £{9) = e'^(Mp) the Souriau mapping ( p| ) associates to 
£{9) the complex number w{9) — e^*^. It follows that we have the identifications 

£{9) = e^''^ and £oo{0) = (e^*^ 26* + 2fc7r), fc e Z. (48) 
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In particular Kp is identified with +1 and (Kp)oo with (l,2fc7r). 

Consider now the tangent plane £{z) at a point z of the Lagrangian subman- 
ifold V. The mapping 

£{.) : V — y A(n) , z ^ i{z) (49) 

is continuous and therefore induces a honiomorphism i^, from the first homotopy 
group of V into that of A(n). In fact, a base point zq ^ V being chosen once 
for all, the mapping 

4 : TTi{V,zo) — > 7ri(A(ri),4) 

(with £o — £{zo)) associates to every loop 7 : [0, 1] ^ F (7(0) = 7(1) = zq) the 
loop ^^(7) of Lagrangian planes defined by ^^(7)(t) = , < t < 1. Using 

the Souriau identification A(n) = W{n) we can associate to the loop ^^(7) in 
A(n) the loop u',7 in W{n) defined by w,,j{t) = w{ii,{j){t)). Lifting that loop 
to Aoo (n) = Woo (n) we get a path 

<i — > (u'»7(t),argdet?«»7(i)) , 0<t<l 

where argdet ?i'*7(t) is a choice of continuously varying argument, uniquely 
determined by a choice of argdet w*7(0). Since u'*7(0) — 1(7,7(1) the quantity 

m{j) = (argdet w*7(l) — argdet w*7(0)) (50) 

ZTT 

must be an integer, only depending on the homotopy class of 7. Formula ( pO| ) 
thus defines a function m : Tri(V,zo) —y Z called Maslov index for loops. The 
integer 771(7) can be intuitively interpreted as follows. Since 7ri(A(n)) = (Z, +), 
A(n) has a "hole". Now, a loop 7 in F induces a loop in A(n), namely the loop 
t t—y £{'yit)) = T^(^t)V, and 771(7) is the number of times ^,7 turns around the 
"hole" in A(n). 

Example 12 Suppose that V is the circle in M^^xMp and 7(t) = g27"t ^ 
< t < 1. We have w^-f{t) = e'*''", < t < 1. The argument of w^-/{t) 
varies from to Att when t goes from to 1; it follows from definition that 
m{'j) = 2. The same argument shows that if j — fJ-i£i + • ■ ■ + fJ-k^k is a loop in 
(S^)'', where £j{t) = e^'^** (Q <t <1) is a loop on the j-th circle, then 

k 

m(7) = 2^A^^.. (51) 

i=i 

The fact that 771(7) is an even integer in the example above is not fortuitous. 
In fact, Souriau Q has proved that: 

V oriented 777,(7) = mod 2 for all 7 £ 7ri(V, Zq) (52) 

(see for an algebraic proof of this property, and the generalization to "g- 
oriented Lagrangian manifolds" ; Dazord Q gives a related cohomological defi- 
nition). 
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Let us next generalize the notion of Maslov index to arbitrary paths in A(n). 
Let ^zqz be a path in V joining zq to a point z and z its homotopy class: z 
is an element of the universal covering V of V . If two paths ^^^^ and ^'^^^ are 
homotopic, then so are their images i*(cizaz) ^"^^ ^*{i'zqz) i'^ ^i^) by 

: F 9 z I — > T,y e A(n). 

This mapping induces a continuous mapping 

iooi-) ■■ V ^ Ao^in) (53) 

which to every z G V with representant "fzoz associates the homotopy class 
£ac{z) of ^,(721,2); obviously the diagram 

V A^(n) 

TT i i-K (54) 

V ^ A(n) 

is commutative (the vertical arrows being the covering projections) . In view of 
the identification A(n) = Win) we can associate to ^^^^ a unique continuous 
path t ^ w{t) (t e [0, 1]) in W{n) such that argdet'u;(t) = a{t), provided that 
we have specified an "initial argument" a(0) for w{0) = i{zo). 

We now impose the following rather restrictive condition on the endpoints 
of the path 72^^: we assume that z is such that 

e{zo) n i{z) = (55) 

and define an "argument function" mg : V — > M by the formula 

1 Ti 

mo{z) = — (a(l) - a(0) + i Tr Log(-w(l)i«(0)-i) + - (56) 
where Tr means "trace of" , and where we define the logarithm by 

Log{~w{l)w{0)-^) = f {[A/ + w(l)(w(0))-i]"'-(A-l)-i/}rfA (57) 

(/ the n X n identity matrix). The right hand side of (|57| ) makes sense in 
view of the following characterization of transversality of Lagrangian planes 
(see (16[ m |6l H): 

Lemma 13 Let £ and £' he two arbitrary Lagrangian planes, and set w = w{£) 
and w' = w{£'). The condition £0 £' = is equivalent to det{w{w')^^ — I) ^ 0, 
that is, to the condition that w{w')~^ has no > eigenvalues. 

We have: 
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Proposition 14 (1) The function mo is integer-valued. It is locally constant 
on its domain of definition {z S F : (-{z) n = O}; (2) mo coincides with the 
function defined in ( p^ when V — and zq — +1; (3) we have for all 7 S 

T^l{V,Zo) 

mo (72) = mo(z) + m(7) (58) 

where m{'~f) is the Maslov index for loops defined by 

Proof. (1) We have, by definition of w: 

exp(TrLog(-w(l)u;(0)"i) = (-1)" det(w(l)w(0)"i) 

= (-1)" (exp(ia(l)) - exp(m(0))) 

and hence exp(27riTOo(i:)) — (— l)"e™^ = 1 so that mo(i) G Z, as claimed. (2) 
If n = 1, y = Si, and zq = 1 then ^(6*) = w{\) = e^*^ and £(0) = w(0) = 1. On 
the other hand the logarithm defined by ( [57| ) is given by in the case n = 1 by 

Lo.9(e-) =^{a~2[^] ^) (59) 

for a ^ Tx (mod 27r), hence 

Log(-z«(l)u;(0)-i) = Log(-e2*^) = i (20 - 2 [^] ^) 

from which follows that mo(i) = [f^/vrj + l, as claimed. (3) Let z be the homotopy 
class of a path 7^^^ and 7 the homotopy class of a loop 7^„^„. Then 7Z is the 
homotopy class of the concatenation 7^^^^ + Iz^z- Formula ( |5^ ) follows, by 
definition ( |56| ) of the Maslov index for loops. 

The last step in the construction of the complete argument index needs the 
properties of the Leray index. 



6.2 The Leray index 

The key to the definition of the Maslov index for paths with endpoints in general 
position is the cohomological index defined by Leray |2^ in the transversal 
case, and generalized by the author jl^ to the non-transversal case. We begin 
by giving a general definition of the Leray index. Recall that Aqo (n) = Woo [n) 
is the universal covering of the Lagrangian Grassmannian A(n). 

Definition 15 A Leray index on {Kaa{n))^ is a mapping 

m : (Aoo(n))^ — > Z 

having the two following properties: (1) the coboundary of m, viewed as a 1- 
cochain, descends to a Sp(n)-invariant cocycle f on A(n); dm = tt*/ (tt the 
projection Aoo(n) > ■^(")J (^J ™ locally constant on each of the subsets 

{(£00,0 :dim(£n£') (60) 
(0<k<n) o/(Aoo(n))'. 
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Condition dm = n* f means that 

mie^,e'^) - m(£^,C) + = fi^J'^n (61) 

and the S'p(n)-invariance of / means that 

f{se, st, si") = f{i, I', e") for all s E Sp{n). 

Notice that the function / automatically is a Z- valued 2-cocycle on A(n): df = 
0, locally constant on each of the sets 

{(£, f, e") : dim{£ D £') = k, dim{£' D £") = k' , dim{£" £) ^ k"} (62) 

(0 < k, k', k" < n). Given a 2-cocycle / on A(n), there exists at most one Leray 
index m satisfying ( |6l| ) (see [l6|). 

We will also need the following simple general property: 

Lemma 16 Suppose m is a real function defined on all the pairs {£ooA'oo) such 
that £0 £' = Q, and such that (|w|) holds for some 2-cocycle f on K{n). Then, 
the formula 

m{£^J'^) = m(^oo, O - ™(C O + /(^, n (63) 

where £'^ is chosen such that £n£' = £n£" defines unambiguously m{£oo, £'oo) 
for all (^00,0 e (Aoo(n))'. 

It is sufhcient to verify that m{£oo,£'oo) independent of the choice of £'^, 
but this is from the cocycle property 9/ = of / (sec |l^, pO[|). 

To every triple {£,£',£") of Lagrangian planes we can associate an integer 
a{£, £',£"), called signature, and defined as being the difference (T-|_ — cr_ between 
the number of > and < eigenvalues of the quadratic form 

Q{z, z', z") = n{z, z') + n{z\ z") + n{z", z) 

on £ © ^' © £" (see |l^, |2^). The signature is an antisymmetric and Sp{n)- 
invariant cocycle: da — 0. Introducing the notation dim{£,£') = dim£ n £' we 
moreover have 

a{£,£',£") = n + ddim{£, £',£") , mod 2. (64) 



Theorem 17 (1) The function m defined by 

m{£oo,£'oo) = ^ (« - + *TrLog(-u,(w')"')) + f (65) 

for £ryo = ('w,a), £'^ = {w',a') with transversal projections: lC\£' = is the 
Leray index associated to the cocycle 

Inert(£, f , £") = ^ {cr{£, £' , £") +n + d dim(£, f , £")) . (66) 
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(Inert is called the "index of inertia" of (i, £' , £")). (2) That Leray index m has 
the following properties: 

m(£oo,0+"^(C^oo) -n + dim(£,0 , TO(^oo,^co) = n (67) 
and the action of TTi{A{n)) on m satisfies 

mix'' ■ £^,X''' ■ C) - "^(^co,0 + k-k' (68) 



where A is the generator of TTi{A{n)) whose natural image in Z is +1 (cf. ^4a))- 
(3) For n — 1 we have, with the notations 

m{e,9')=[^]+l. (69) 

Proof. We first remark that Inert (^, £', £") always is an integer in view of (Q). 
We have shown in |jll|, ^ that the function ^ defined on all {(f oo, ^'oo) ■ £n£' = 0} 
by /Lt = 2m — n (m defined by (|65|)) satisfies 

iii£oc,iL) - A^(^oo,c) = ^{i,i\n- 

It follows that 

mi£oo, O - ™(^oo, C) + mi^LJ'L) = \ H^, n + n) 

if the planes £, £' , £" are pairwise transverse. Since in this case d d\in.{£ ,£',£" ) = 
0, the existence of m follows from Lemm a since Inert obviously is a Sp{n)- 
invariant cocycle. Formulas (^, (|6^), (|68D are obvious consequences of (^, 
( |65| ) when £ n £' = 0, and of ( ]6l| ) in the general case since £oo and k ■ £oo have 
same projection £ on Lag(n). Let us finally prove property (3). Suppose first 
that — 9' (mod tt). Then ( |6^ ) immediately follows from (|65|). Suppose 
next that - 9' = kn. Choosing ^ such that £{0) n ^(6*") = 0, we have 

m{9, 6') = k + Incrt(^(6l), £{9), £{9") = k+l 

in view of (|63[), concluding the proof. 



Remark. There is a deep and interesting connection between the Leray index 
m and the Maslov index on the metaplectic group Mpin) (i.e. the unitary 
representation of the double cover Sp2{n) of Sp{n)); see 

6.3 Definition of the Argument Index 

We now have developed the machinery we need to define the complete argument 
index generalizing (|2l|)~(p^). 

Consider again the mapping £oo(-) : V — > Aoo(«) defined by (p4[). We 
denote by £a,oo an element of Aoo(n-) with projection £a G A(n). 
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Proposition 18 The function nia : V — > Z defined by 

nia{z) = m{£oo(z),£a,ryo) (70) 

has the following properties: (1) Suppose that = io and ia.oo ~ ^o,oo is the 
homotopy class of the constant loop with origin Iq- Then ma{z) is given by ( [5^ 
if£{z)n£o = 0; (2) We have 

rriaijz) = rriaijz) + 771(7) (71) 

for all 7 G T^iiV) and z G V ; (3) We have ma{z) — m{9) when V = , 
£a — and zq ~ +1. 

Property (1) is obvious in view of ( p6|) and ([6^). Property (2) follows from 
property ( |68|) of the Leray index. Property (3) follows from part (3) of Theorem 

The following result makes explicit the effect of a change in base point: 

Proposition 19 Let nia, mp be the Maslov indices associated by ( |5^ to arbi- 
trary elements £ a ^00 o-nd £p ,^ of K^{n). We have 

ma{z) - mp{z) = m{£a^oo, £p,oo) - Inert(^a, £p, £{z)) (72) 

where z is the projection on V of z Cz V . 

Formula ( j?^ ) is of course an immediate consequence of property (|6^) with 
the choice / — Inert. 

We next define the waveforms on Lagrangian manifolds. 

7 Waveforms 

We set out to generalize the construction of the square root of a de Rham 
form on the circle, as outlined above, to the general case of a n-dimensional 
Lagrangian submanifold V (which we again suppose connected, but not neces- 
sarily orientable). We will use the following standard notation and terminology: 
the caustic of V is the closed set 

C ^ {zeV ■.£{z)C^W^^Q) . 

More generally, we will call "caustic of V relatively to the direction £a G A(n)" 
the closed set 

Cc.^{zeV : £{z) r\£c,^0} 
and we denote by its complement V\Ca'- 

Va^{zeV : £{z) n£a = 0} . 
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7.1 De Rham Forms and their Square Roots 

Remark. The introduction of de Rham forms in SM should not be too sur- 
prising after all. It is well-known in Physics that many phenomena exhibit this 
dependence on orientation, the most elementary example of this phenomenon 
being the magnetic field, which is a "pseudo-vector" (see the lucid discussion in 
Frankel's book |^). On the other hand, the necessity of the inclusion of half- 
densities (or half -forms) in quantum mechanics has been remarked a long time 
ago (it is of course immediately suggested by Van Vleck's formula Q (also see 
[|8[ ^). (Historically, the systematic use of these objects goes back to the work 
of Blattner, Kostant and Sternberg (see ||, |2^, |39|]).) 

A m-density p G \n"^{V)\ on V (to G M) is a smooth section of the line- 
bundle |A™(rF)| of 1-densities on TV. Recall that by definition every p{z) G 
\A"'{TzV)\ is a mapping 

p{z) : T^V X • • • X T,y — > C 



such that 



p{z){ui, ...,Un) = |det p{z){Aui, ...,Aun) 



for every A G GL{n,C) and all vectors iti,...,u„ in TzV. Let {Ua,fa)a be an 
atlas of V. The local expression of p G |0"'(y)| in each chart {Ua,fa) is 
p^{x) \dx\"^, the functions p^ G C°°{fa{Ua)) satisfy the matching conditions 



, . dfaf3 , ^ 

det^(.) 



Pp{x) ,xef0{Uc.nUp) (73) 



where we have set f^p — fa° fp^ ■ In particular, if and are two coordinate 
systems on [7^, then 



p^{x) \dx\' 



d t 

dx' 



Pa{x') \dX 



I I ^ 



if a; = fa{z), x' = /^(z). 

Suppose now that to = 1; we write |Sl^(y)| = |f^(V)| and call elements 
of |r2(y)| simply densities. Obviously, each p{z) G |A^(r2l^)| is homogeneous 
with respect to scalar multiplication, but it is not additive. The definition of 
de Rham forms reinstates additivity: a de Rham form p, G ^r{V) associated to 
a density p G |f^(V)| is a smooth section of the line bundle h.{V) obtained by 
assigning to each p{z) G |A(r2F)| the mapping 

p{z) : X • • • X — > C 



defined by /i(z)(ui, ...,u„) = if the vectors wi, ...,Un are linearly dependent, 
and by 



p{z){B^{z))^±p{zmz)) 
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if they form a basis B{z) of £{z) = TzV; the notation ± refers to whether this 
basis is positively or negatively oriented relatively to the orientation at z defined 
by a local chart {U, /) at z. Due to the inclusion of the factor ±1 in its definition, 
niz) is linear, and antisymmetric. 

Let us now be more specific, and assume again that y is a Lagrangian 
submanifold. We denote by fa the orthogonal projection of Va on ^q,; it is a 
local diffeomorphism onto its image, so that each orientation of ia determines 
an orientation at z d Va- Let now U be an open neighborhood of z in V^. 
Choosing U sufficiently small, the pair {U, fa) is a local chart at z. The open 
set U is orientable, and each of its orientations is determined by the choice of 
an orientation of £a, that is by the datum of an element ^Q,of the double cover 
A2(n) with projection £a- The restriction £u{') of the mapping z £{z) to U 
can be lifted to two continuous mappings z i—> £f}{z) G A2(n), corresponding 
to a continuous positive (resp. negative) choice of orientations of the tangent 
planes. Let p E \ft\ (U) be a density on U, and (-ij ra^z) be two elements of 
Aoo(n) with projections £^{z) G A2(n). Let £a,oc be an element of Aoo(n) with 
projection £a on A2{n). We claim that: 

Lemma 20 The formula 

^^(z)(S±(z)) = (74) 

defines a de Rham form on U . 

In fact, if we change S+(z) into B^{z), then we have to change ^t/oo(-^) i^^o 
£(}oci^) ™ formula (|7^). Since both ^^oo(^) ooi^) have same projection 

£{z) g A(ri), we must have £^ ^{z) — X'^ ' ^jj oo(^) ^'^^ some integer k and hence. 



by (68) 



"^(4,00 ^a,oo) = m{£^^^{z), £a.oo) + k . 

Now, the integer k must be odd, because if it where even, then £'lj ^{z) and 
ooi^) would have same projection £^{z) on A2(n). Thus 

^Jiu{z){B+{z))^-^Jiv{z){B-{z) . 

Similarly, if we reverse the orientation at z, that is, if we replace £a by an 
element of A2(n) defining the reverse orientation, then we must replace £a.oo by 
A • £a.,oo, which leads to replace m{£f} ^{z), £a,Qc) by 

and thus again reverses the sign of piu{z){B^{z)). The lemma follows, noting 
that the mappings z i—> £f} ^{z) arc locally constant, and hence smooth. 

Formula ( [z^ allows to define locally the argument of a de Rham form by 

HTgHjjiz) = ™(4„^(z),£+^)7r (75) 
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and hence the square root of by the formula 

^{z){B^{z)) = (^^■'°-Vp(^)(^(^)) (76) 

It turns out that this formula can not generally be extended to define a 
global argument for a de Rham form (cf. for example the density r\dO\ on the 
circle S}.). However, we are going to show that this can always be done outside 
the caustic Ca relative to provided that we work on the universal covering 

of y. 

Proposition 21 Let ^ he a de Rham form on V , associated to a density p £ 
For every Va there exists a choice of ia.oo G ^ooin) such that the 
restriction of ^ to Va is given, for z g T^~^{Va), by 

= (77) 

where m'^{z) = m{ioo{z),£a,co) and m~{z) = m{iooiz), X ■ ia.oo)- We can thus 
define the square root of p, of Va by the formula 

V7^(i)(S±(z)) ^ t^^^'^^K^WW)- (78) 

If z e Fq n Vfj, then 

yM^(z)=z"°^WVM^(z) (79) 

where the function ra^p :VaC\Vp^'Lis given by 

7nap{z) = m(ta,oo, ip^oo) " Inert(£Q, i{z)). (80) 

Remark. We have introduced similar notions in |l^ ; however the choice 
Ma = i"^"^'^^ ^/p for the square root of a half-density was postulated in a rather 
ad hoc manner. 

It is instructive to interpret the constructions above in terms of the oriented 
double covering V of the manifold V. Recall (see for instance X, §4) that 
V is constructed, for an arbitrary submanifold V , in the following way: let 
(Ua, fa)a be an atlas, and define, for Ua^Up 0, locally constant mappings 
gap -UanUp — > {+1,-1} by 

gap{z) = DfafiiffAz)) \Dfap{fp{z))\-' . (81) 

where the fap = fafp^ are the transition functions. Evidently gapgp^g-ya{z) = 
z for z £ Ua C] Up n , hence there exists a twofold covering n : V — > V with 
trivializations 

Ta ■■ n-\Ua) -^UaX { + 1,-1}. 

The transition functions 

Tap = TaT^^ : Up X { + 1, -1} > [/„ X { + 1,-1} 
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are given by rQ^(z, e) = (z, (/^^(z)) for z e UafMIp and e = ±1. This allows one 
to construct an orientable atlas {Ua.e, fa.e)a,e of V by setting Ua^e — T~^{Ua x 
{s}) and defining fa^e '■ Va,s — > by the formulas 

r Uz) for e = +l 
/„,,(T-i(z,e))= <^ 

[ <yfa{z) for £ = -1 

where cr is M" — > M" changes, for instance, the coordinate xi into —xi and 
leaves the other coordinates unchanged. One has the following property: V is 
orientable if and only if the double covering V is trivial: V = V x {+1,-1}, 
and V is connected if and only V not orientable. 

When V is Lagrangian, we have the following interesting result that shows 
that the oriented double cover can always be identified with the product V x 
{+1, —1} (but of course not equipped with the product topology when V is non 
orientable!): 

Proposition 22 Suppose that V is not orientable. Then, each of the mappings 

^a.:V X {+1,-1} 

defined by = (z, wh ere z (z V has projection z Cz V , is a 

bijection. The restriction ^aof to Va = {z : i{z) D ia} is a homeomorphism, 
and the transitions — ^a^^^ o.i'e given by 

^ap{z) = (z, to(£„,oo, ^/3,oo) - Inert(4, Ip, t{z)) (82) 

for 5 e t4 n %. 

Proof. We first note that $0, is well-defined: if z' and z both have projection 
z e V", then z' — 7Z for 7 e ■ki{V), and hence 

ma{z') = rUai^z) = mct(z) mod 2 

in view of (|5^) since V is orientable. A similar argument shows that each 
mapping $q, is injective: if $0(2') = $q(z), then z' = z and z' = 7Z with m^j) 
even, so that 7 G Tri{V), and z' = z. To prove that $q is surjective, it suffices to 
note that if $q-(z) = (z,e), then <i>o,(z') = (z,— e) where z' is the projection on 
V of 7z, z has projection z, and 7 G tti{V) is such that 771(7) — 1 (the existence 
of such a 7 follows from (^2|) since we are assuming V non-orientable). Finally, 
4>Q, is locally constant on V^, and hence continuous; formula (|2|) follows from 

&■ ■ 

7.2 Definition of Waveforms 

Let if be the phase of the Lagrangian submanifold V, and /i a de Rham form 
associated to a density p on P^. 
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Definition 23 A waveform on V is the datum, for each £a G of an ex- 

pression 

^^{z)^eTi'^^'^^iz) 
where the /z^ are associated to a same density p on V . Equivalently, 

Defining the action of 'ni{V) on by ^'^{z) — ^'(72), we say that is 
defined on if 7^' = for all 7 € t:i{V). 

The following result relates our constructions to the minimum symplectic 
area postulate: 

Proposition 24 (1) A waveform is defined on V if and only if V satisfies the 
Maslov quantization condition 

If 1 

2^ J pdx+ -m.(7) is an integer (83) 



for every loop j in V . (2) When V is oriented, the condition ( pq j reduces to 
the condition 

^JpdxeZforalljeTTiiV). (84) 



Proof. The second statement of the theorem follows from the first in view of 
(p2|). By definition of a waveform we have to prove that the condition 7^'(i) = 
\l/(7z) is equivalent to (^3|). In view of (|2|) and ( [7l|) we have 



*(z) 



7^1' (i) = exp 
hence 74* = ^ is equivalent to 

If TT 

— / pdx H — ^(7) = mod 27r 
h J J 2 

which is of course the same thing as ^ 

7.3 The Hamiltonian motion of waveforms 

The waveforms we have defined are time- independent, they are thus adequate 
for the study of stationary processes. However, if we want to use them for the 
dynamical study of quantum systems we have to define how Hamiltonian flows 
act on them. 

Consider an arbitrary function H — H{x,p,t) defined on some open subset 
D X Rj of the extended phase space R" x x Mt; we denote the time-dependent 
flow of Xh - {VpH, ~V^H) by ft.t'- 
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We need the following property of the Leray index. Let Spoo{n) be the 
universal covering of Sp{n); the usual action Sp{n) x A(n) — > A(n) induces an 
action Spoa{n) x h-oa{n) — > Aoo(n). We claim that following essential formula 
holds: 

to(Soo^oo,SooO = "^(^00,0 (85) 

for all (sooj^oo) G Spoo{n) x Aoo(f^)- 

Formula follows from the fact that the cocycle / associated to a Leray 
index is S'p(n)-invariant and from the uniqueness of the Leray index associated 
to a given cocycle (see Lemma |l6|): both mappings (€00, ^'oo) 1 — m(£oo,^'oo) 
and (^ooj^'oo) ' — "^(scxj^oo, Soo^'cx)) satisfy the conditions in definition |l5|, and 
are hence identical. 

The Jacobian matrix st^t'{z) of ft^t' at z being symplectic, we can lift the 
mapping t ^ st^t'{z) € Spin) to a mapping 

t ^ istM{z))oo G Spoo{n) 

such that {st,t(,z))oo is the identity of Spoo{n)- Setting 

■mo{z,t) = m(^o,oo, (st,t'(2;))oo-^(2)) (86) 

we then define the value of "if at time t by the formula: 

^-(2,0 = e^'^(*^*)z'"«(*'*)yC/^(z). (87) 

Let /t^j/ be the mapping which to ^'(i,t') associates ^(i,t). These mappings 
satisfy the Chapman-Kolmogorov condition: 

ft.t' ft',t" = ft.t" ■ (88) 

as immediately follows from the fact that 

{ft,t'ft',t")*P ~ ift,t')*ift\t")*P- 

7.4 Shadows, and their Motion 

Suppose that the Lagrangian submanifold y is a simply connected graph p = 
S/x^ix). V is then automatically quantized, because 771(7) — for every loop 7 
in V. We denote by S the operator which to every waveform "ifonV associates 
the coefficient of its local expression in the chart {Vt , ttx ) where ttx is the 
projection {x,p) 1 — > x on configuration space. Thus, if ^' has local expression 
e^*(^)a(x)|da;|i/2 in {V.iix), then 

5(*)(x) = ei*(")a(x) . 

We will call S(5') the shadow of . Suppose now that V is an arbitrary La- 
grangian submanifold (i.e. we relax the conditions that be a graph, or simply 
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connected). We moreover assume that the quantization condition (p3|), holds 
for V. Then, given a point x there wiU in general be several charts {Uj,iTx) 
such that X e TTx{Uj). In this case we define the shadow of a waveform ^ as 
being the sum 

3 

calculated at the point (x, x'j, t, t'), and rrij is the Morse index of the trajectory 
from x'j to x: it is the number of conjugate points along that trajectory, that 
is, the number of points where detHess2:_2^'(5) is zero, or infinite. 

Write now the wavefunction at time t' in the familiar "oscillatory" form 

where $ and a are smooth functions, a > 0, both defined for {x,t') e X x Mt , 
X some open subset of M" (we do not assume that X is simply connected. For 
fixed t' the function ip{-,t') is the local expression of a Lagrangian waveform 
^{.,t') on the graph Vf of the gradient of the phase <i>(.,t'): 

Vt>^{ix,p):p = V,M^,t')}. 

In fact, 

^^>{z,t') = et"^**'*'^ /* (a(x,t') |d"a;|'/') (89) 

where /* (^a{x) \d"x\^^'^^ is the pull-back to Vf of the half-density a{x) 
on R^. 

Let us finally relate our waveforms to the approximate solutions to Schrodinger's 
equation studied by Maslov and Maslov and Fedoriuk . Writing the ini- 
tial wavefunction as 'i'o{x) = exp [■|S'o(a;)], these authors propose expressions 
of the type 

i 

where x'j, Sj and fij{x,t,t') are defined as follows: let again {ft.t') be the time- 
dependent fiow of H, and denote by V* the image by ft.t' the image of the 
Lagrangian submanifold Vt> : p = VxS{x,t'). Given a point x of the projection 
of Vt on M" there is (under adequate assumptions on U and ) a finite number 
of points {x,pj) e Vt and {Xj,Pj) G Vf such that {x,pj) — ft,t'{x'pp'j). The 
phase Sj{x,t) in ( pO| ) is then given by the integral 

Sj{x,t) = S{x,x'.i;t,t') = / pdx — H ds (91) 



dx 



-i/z 



exp[f5,(x,t)] -^oixj^t') 



(90) 
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calculated along the trajectory leading from x'j at time t' to x at time t. The 
integers iJ,j{x,t,t') in ( |90| ) are the Morse indices of these trajectories; these 
indices are obtained by counting the number of conjugate points along each 
trajectory. It turns out that for short time intervals t — t' formula ( |90| ) reduces 
to 

-1/2 

(92) 

where S is just the classical action function evaluated from {x',t') to {x,t). In 
fact, if t — t' is sufficiently small then Vt will be a graph and there will exist 
exactly one point x' such that {x,p) = ft^t'ix' ,p') for {x,p) € Vt- (Formula ( p2[ ) 
was actually already discovered in 1928 by Van Vleck p^.) 

Conclusion. We have achieved our goal, which was to construct a scmiclas- 
sical mechanics based on a topological principle without any reference to the 
usual semiclassical approximations based on the WKB method (from which SM 
historically originates). Semiclassical mechanics thus appears to be a theory in 
its own right. We have not examined in which sense our theory approximates 
CM or QM, nor have we given any applications. As is the case for CM or 
QM, the domain of validity of SM can be determined by experience. It would 
certainly be interesting to develop further the consequences of the minimum 
symplectic area postulate in the following directions: 

(1) The quantization of non-intcgrablc systems; it is well-known that periodic 
orbits play a fundamental role in such systems (see e.g. ^); the theory 
of Lagrangian paths as outlined in Section ^ and further developed in pof is 
certainly useful in this context (this theory gives a mathematically rigorous 
justification of the recent constructions in Sugita pst ; 

(2) Statistical mechanics and thermodynamics: the minimum symplectic 
area postulate could be used to push further phase space "cell quantization" as 
we outlined in [|l9j . Formula (^ relating symplectic area and volume shows that 
the volume of such cells corresponding to a quantized ball B{\/H) is /i"/2"'n!, 
which is consistent with Bose-Einstein statistics. 
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